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I calculate the superfluid density of a non-equilibrium steady state condensate of particles with 
finite lifetime. Despite the absence of a simple Landau critical velocity, a superfluid response sur- 
vives, but dissipation reduces the superfluid fraction. I also suggest an idea for how the superfluid 
density of an example of such a system, i.e. microcavity polaritons, might be measured. 
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The observation of superfluidity is one of the most com- 
pelling signatures of quantum coherence in systems such 
as Helium and cold atomic gases [if. Recently, there 
has been much activity exploring condensation of mixed 
matter-light excitations, i.e. semiconductor microcavity 
polaritons (see Q for a review) as well as recent exper- 
iments on photons in dye filled cavities Microcavity 
polaritons consist of superpositions of excitons confined 
in quantum wells, and photons confined in semiconduc- 
tor microcavities. Because of the photonic component of 
these particles, they have a finite lifetime, and so any con- 
densate will be a non-equilibrium steady state, where loss 
is balanced by injection of new particles. As well as these 
matter-light systems, recent experiments on continuous 
loading of cold atoms into traps [H suggest that simi- 
lar questions of superfluid properties in non-equilibrium 
steady states may soon be accessible in cold atoms. 

Finite particle lifetime prompts questions about the 
meaning of superfluidity when the particles involved in 
the condensate are continually being replaced. To see 
why such questions arise, one may observe that finite 
particle lifetime changes the excitation spectrum from 
the linearly dispersing Bogoliubov sound mode to a dif- 
fusive mode 0, \ 



tices 111 have been observed. While these various ex- 
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of the spectrum of long wavelength modes is commonly 
invoked in explaining the existence of superfluidity, i.e. if 
the low energy excitations are linearly dispersing Bogoli- 
ubov sound modes, then there exists a critical velocity, 
such that for a fluid flowing below this velocity, it is not 
energetically favorable to create quasiparticle excitations, 
and thus the flowing superfluid remains stable [l[ . If one 
naively defines a critical velocity by the dispersion of the 
real part of w(fc), one finds that for the non-equilibrium 
spectrum, such a critical velocity vanishes. As has been 
pointed out elsewhere @, Sj, a more careful argument 
regarding the response to a static defect indicates that 
there may still be a particular velocity at which there 
is a sharp onset of drag. Nonetheless, the above illus- 
trates why non-equilibrium condensates require one to 
re-examine questions of superfluidity. 

Experimentally, aspects of superfluid behavior have 
been studied in microcavity polariton systems, and fea- 
tures such as quantized vortices 0, suppression of scat- 
tering off disorder [l(| and metastability of induced vor- 



pcrimcnts show that aspects of superfluid behavior can 
be seen in such non-equilibrium condensates, they leave 
open an important question, namely what is the super- 
fluid fraction. Even in superfluid Helium, at non-zero 
temperatures there are drag forces due to the normal fluid 
component. However, the normal and superfluid compo- 
nents can be clearly distinguished by their response to a 
slow rotation [l|: At low angular velocities, the superfluid 
cannot rotate, so only the normal component rotates, 
thus reducing the effective moment of inertia. Determin- 
ing the normal and superfluid densities thus gives a fuller 
description of the superfluid properties of a system than 
does a binary distinction between superfluid and non- 
superfluid systems. As such, the aim of this letter is to 
discuss the normal and superfluid densities of an open 
dissipative condensate, and to suggest how this might be 
explored in the microcavity polariton systems. 

The same superfluid density as introduced above can 
be found from the current-current response function 
Xij(q)- This response function gives the particle current 
Ji{o) = J2k iPk+qliCZk + q)ipk (where the current vertex 
7i(fc) = hi /2m and tl = 1 throughout) due to a pertur- 
bation h h -J2 q fi(q)Ji(q), Le - Mq) = Xtj(q)fj(<j 

In an isotropic system, one may then define p s , p n as 



PnO~i. 



(i) 



The superfluid component picks out and responds only 
to the irrotational (non-transverse) part of the applied 
force. The advantage of this approach is that it allows 
an explicit calculation for the open dissipative system, us- 
ing the Schwinger-Keldysh approach for non-equilibrium 
systems (see e.g. fl2|h 

In the following, I will consider the simplest model of 
a weakly interacting dilute Bose gas: 
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where = k 2 /2m. In addition, one must include pump 
and decay processes such that the bare inverse retarded 
Green's function = uj — ej, + in — ip{uS) where 

the pump has the form p(uj) =7 — 770; and k describes 
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FIG. 1. Types of Feynman diagram required for the re- 
sponse function to one-loop order. Straight lines indicate 
non-condensate excitations. Filled symbols involve the con- 
densate, arising from either interactions (circles) or coupling 
to currents (squares). Wavy lines indicate source fields cou- 
pling to the current vertices. 



decay [13j- This form of pumping is motivated by recent 
works b y W outers and Carusotto Q (as well as related 
models [14j), and simplifies the calculation of Xij com- 
pared to models with density-dependent pump processes. 
When condensed, such a model has the diffusive spec- 
trum £(fc) = — ier + \/c 2 fc 2 — a 2 with a = T]fj,/(1 + ?7 2 ), c = 
\J /i/[m(l + r/ 2 )]. Thus, finding a non-zero supcrfluid 
density in such a model addresses how the diffusive spec- 
trum affects superfluidity. 

To correctly find the supcrfluid response function 15 1 
requires vertex corrections in the response function. In 
equilibrium, this can be avoided by using sum rules that 
result from conservation of density, however with finite 
particle lifetimes, this is not necessarily a-priori justified. 
I will postpone until later the discussion of how the vertex 
corrections are to be determined and next summarize the 
physical reason that a superfluid density can survive. 

Fig. [1] illustrates the classes of Feynman diagram (in- 
cluding vertex corrections) that result at one loop order. 
The first five diagrams contribute to the superfluid den- 
sity, while the last gives the normal density. This can 
be seen by noting that the first five diagrams all have 
the current vertex scatter a particle out of the conden- 
sate, and thus involve a factor ji(q) cx qi, hence they all 
contribute to Xij lilj- I n order that the supcrfluid 
density does not vanish, it is crucial that the fluctua- 
tion propagator D R (q,uj = 0) that also appears in these 
five diagrams behave as 1/q 2 at q — > so that overall 
Xij oc qiqj/q 2 remains finite. The existence of supcrfluid 
density therefore depends on how the denominator of the 
Green's function behaves. 

In thermal equilibrium, the Green's function behaves 
as D R (q,ui) cx [(w + iO) 2 — (?q 2 \~ x and so the correct 
scaling of D R (q,uj = 0) is dependent on having the a 
linear spectrum, hence the relation of the Landau crit- 
ical velocity and supcrfluid density. However, despite 
the changed spectrum of the open dissipative system, 
one has D R (q,u>) cx [uj 2 + 2i<roj — c 2 g 2 ] -1 and so the 
Green's function at u = still scales as D R cx 1/q 2 , 
yielding a non-vanishing supcrfluid density. Such behav- 
ior of the Green's function has also been seen to exist 
in several other models of non-equilibrium polariton con- 
densates j^,|f|. The fact that this structure of the Green's 
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FIG. 2. Normal density vs temperature for a variety of pump 
and decay rates and obeying the equilibrium fluctuation dis- 
sipation theorem. Relative pump and decay rates are scaled 
so all three curves correspond to the same condensate density, 
while the influence of pump and decay varies. 



function leads to a supcrfluid density, despite the modi- 
fied spectrum, is the first main result of this letter. 

A second result is the effect of finite particle lifetime 
on the normal density. In an equilibrium single compo- 
nent system, the normal density vanishes at zero tem- 
perature [l| . The normal density of the non-equilibrium 
system can be straightforwardly calculated since, just as 
in the thermal equilibrium case, there are no vertex cor- 
rections at one loop order [15j, so one finds (in 2D): 
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where the Green's functions and Pauli matrices cr, are 
written in Nambu space, i.e. D R (t,t') = —i6(t — 

0([* fc (*) > *fc(*')]> > *I = (fk^-k)- Even for a ther- 
malised case, using the equilibrium fluctuation dissipa- 
tion theorem, Df (uj) = (2n B (uj) + l)[D R {u) - D£(w)\, 
one finds that the presence of pump and decay terms af- 
fect the normal density. As shown in Fig. [2l the normal 
density does not vanish at zero temperature. 

Having shown that supcrfluid density need not vanish 
in a dissipative condensate, but is reduced by finite life- 
time, one may then ask how the superfluid and normal 
densities could be measured in such a system. As an 
illustration, the following suggests a method that uses 
the polariton polarization degree of freedom [lrjj in or- 
der to apply ideas that have only recently been proposed 
for how one might measure of supcrfluid density in cold 



atom systems 17|, LL8(- A number of alternative meth- 
ods likely also exist, such as adaptations of proposals to 
create gauge fields in coupled photonic cavities [Uj]. To 
adapt the approach in [18|, one first considers the effec- 
tive Hamiltonian Hb due to an inhomogeneous real mag- 
netic field, in the space of polariton polarization states. 
If the splitting of polarization states is always large, one 
may restrict to the adiabatic ground state Because 
the polarization composition of this ground state varies 
in space, there can be a non-trivial synthetic gauge field 
in this subspace gj4 synt h = i('0|VV')- Thus, a real mag- 




FIG. 3. Panel (a), microcavity sample (/xc) placed between 
(imbalanced) anti-Helmholtz coils to induce a field B roa i. 
Panel (b), velocity and angular velocity vs radius for A Bynt ^ 
as discussed in the text. 



nctic field acting on polarization degrees of freedom can 
induce an artificial vector potential acting on the (neu- 
tral) polaritons. This artificial gauge field can mimic a 
rotating frame, thereby allowing one to distinguish the 
superfluid and normal response to rotation 

In order to illustrate how this might work, one may 
consider the real magnetic field produced by an im- 
balanced anti-Helmholtz configuration, as illustrated in 
Fig. Eta), so that the magnetic field in the microcav- 
ity has the form B = ((3x, f3y, B z ) for small in-plane 
coordinates x,y with constant B z . Microcavity polari- 
tons typically involve heavy-hole excitons, so that po- 
lariton polarizations ±1 imply electron and hole spins 
(=Fl/2, ±3/2). As such, while a B field perpendicular to 
the microcavity simply splits these polarization states, 
the in-plane field is more complicated, as it mixes the 
polariton states with non-radiative excitons with spins 
(±1/2, ±3/2). Furthermore, depending on the crystal 
symmetry and quantum well growth direction, the lead- 
ing order coupling between ±3/2 hole states may either 
be linear or cubic in magnetic field fiol ] . In order to illus- 
trate the basic idea, I will avoid these complications, and 
consider the simplest situation, where a linear coupling 
exists [2lJ]. Adiabatically eliminating the non-radiative 
excitons, the effective Hamiltonian for the polariton po- 
larization is H = X[£ 2 a z + r 2 (e 2i,#, cr_ + e~ 24< V + )] where 
re 1 ' 1 ' = x+iy and the length i encodes the ratio of in-plane 
and perpendicular fields. One may then find the ground 
state \ip), ensuring this is smooth as r — > 0, and thus find 
qA = i{ip\Vtp) = (<j>/r)(l-l 2 /Vr 4 + l 4 ). This gauge field 
is equivalent to a rotating frame with qA synt h = mu x r 
hence Lu(r) = qA^/rar. This is shown in Fig. |3jb). 

Detection of the response to this rotating frame could 
potentially be done by imaging the momentum and en- 
ergy distribution of the polaritons [i| . If the condensate 
is concentrated around r ~ t, one has the maximum ve- 
locity, corresponding to an energy shift for the normal 
component of Ai? max = (l/2)mw^ ax = 0.08/to^ 2 . For 
i ~ 0.5/im, and m = 10 ro e i ec tron this corresponds to 



0.2mcV. Observing the differential shift of luminescence 
as one varies £ by varying B z could allow one to extract 
the superfluid fraction. 

I now turn to discuss in more detail how the super- 
fluid density can be calculated in the Schwinger-Keldysh 
approach, and to explain the origin of the vertex cor- 
rections in Fig. Q] and the form of the normal density 
in Eq. ([3]). Following the path integral approach to the 
Schwinger-Keldysh formalism 12J , the response function 
can be written in terms of a generating functional as: 



d 2 Z[f,9] 



2 dfMdOji-q) 



V(4>^)e l(s+5S) . (4) 



Here, S is the Keldysh action |12| due to the Hamil 



tonian in Eq. ([2]) and the pump and decay terms S 
E ^fc^Vfc' - II dtd 2 r%$ q $ cl W 2 cl + tf) + H.c] where 
the fields ip are written in the Keldysh space for clas- 
sical/quantum fields [l2|, and the retarded, advanced 
and Keldysh components of the inverse Green's func- 
tion Vq 1 arc as discussed above. The source term 
SS = Y J '^k+q[Tifi{q) + {Tz J riT2)Qi{q)\ r )i{2k+q)'4) k , where 
Pauli matrices Tt are in the Keldysh space. The presence 
of two separate fields is necessary due to need to calcu- 
late a normal ordered current (derivative with respect to 
9) linearly dependent on a classical force (field /). 

In the non-condensed state, one may immediately de- 
termine Z at leading order by neglecting interactions, 
and performing the (Gaussian) integration over fields 
(4>,i>)cLq- This yields m\ij = pnkj, with p n given by 
Eq. ©, where the Nambu structure in the normal state 
is trivial. When condensed, vertex corrections become 
important. These vertex corrections can be found by 
use of an "honest saddle point" of the partition func- 
tion, i.e. determine the saddle point in the presence of 
the source terms /, 9, and then integrate out fluctua- 
tions about this new saddle point. One thus finds a form 
Z cx exp{iS* [/,6>] - ±Trln(l +VA[f,6})} where S is 
the saddle point action, T> are the Green's functions (for 
/ = 6 = 0) and A is the self energy due to the fields 
/, 9. Both So and A involve terms arising from the shift 
of the saddle point field ip in the presence of the source 
terms /, 9, and thus both So and A have higher order 
contributions of /, 9. One may then expand these terms 
to quadratic order in /, 9 and evaluate \ij y i & Eq. (HI) . 

Taking derivatives with respect to /, 9 (indicated by 
primes and subscripts), one finds Xij = 5^0 / 4 e- + 
iTr(VA'f. 9 .) - iTt(VA' f VA' e .). Comparing the three 
terms in this expression to the diagrams in Fig. [T] the 
first diagram arises from the first term, the second two 
diagrams arise from the second term, and the last three 
diagrams from the third term. After explicitly evaluating 
the shifts to the saddle point, and the self energy A, one 
finds the explicit form: 
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where iV Qb - / ^Tr(D§ +q a a D^ a b 
terms of Pauli matrices cr^ in the Nambu space and A 
C/|V'o| 2 - The terms in Eq. ((5J are arranged in the same 
order as the corresponding diagrams in Fig. [TJ 

A number of technical issues regarding regularization 
are worth noting. Firstly, as is known elsewhere (see 



e -g- 123 an d refs. therein), there can be cases where it 



is necessary to return to the discrete time coherent state 
path integral in order to correctly incorporate causality 
in performing integrals. The current problem is such a 
case, and the term arising from this is the 1 on the first 
line of Eq. ([5]). Secondly, in order to give an ultravio- 
let finite expression, it is necessary to perform the stan- 
dard T-matrix regularization of the contact interaction: 
U- 1 -> U~g -mj f^[2(e + /i)]" 1 . Thirdly, the appar- 
ently singular terms involving qiqj/q 2 e q in fact cancel, 
leaving only finite contributions as q — > 0. One may ver- 
ify that Eq. ([5]) recovers the expected equilibrium result 
in the absence of pumping and decay. 

One may note that the expression in Eq. ([5]) does not 
explicitly involve details of the pumping. This is because 
the only nonlinearity included being the interaction term 
U. This means that Eq. (|5|) survives for general models 
of pumping and decay, and so is more generic than the 
particular model of pumping used to derive it. 

In conclusion, the supcrfluid density of a non- 
equilibrium open dissipative condensate need not vanish, 
despite the non-existence of a Landau critical velocity. 
This is because the poles of the response function, which 
give the spectrum, do not uniquely determine the form 
of the response function at zero frequency, which is the 
quantity that defines the superfluid density. Such a su- 
pcrfluid density could potentially be measured in a po- 
lariton system by using real magnetic fields to engineer 
an effective rotating frame. The current-current response 
function can be explicitly calculated using the "honest 
saddle point" approach. Such an approach would also 
allow calculation of dynamical response functions, allow- 
ing a more nuanccd understanding of the distinctions be- 
tween static and dynamic supcrfluid phenomena in open 
dissipative condensates. 
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